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We calculate pseudorapidity (rj) asymmetry in pA and dA collisions in a pQCD-improved parton 
model. With the calculations tuned to describe existing spectra from pp collisions and asymmetric 
systems at midrapidity and large rapidities at FNAL and RHIC energies, we investigate the roles of 
nuclear shadowing and multiple scattering on the observed asymmetry. Using this framework, we 
make predictions for pseudorapidity asymmetries at high pr and large rj in a wide range of energies 
up to LHC. 
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I. INTRODUCTION 

Collisions of asymmetric nuclear systems, like proton- 
nucleus or deuteron- nucleus collisions, attract significant 
experimental and theoretical attention at present. At the 
Relativistic Heavy Ion Collider (RHIC), Run 8 included 
deuteron-gold collisions at 200 AGcV with a luminosity 
increase of about an order of magnitude compared to the 
deuteron-gold run of 2003. One of the physics goals was 
to provide a high-statistics "cold" nuclear matter data 
set to establish a definitive baseline for "hot" nuclear 
matter, created in gold-gold collisions. The benchmark 
role of the "the deuteron-gold control experiment" for 
e.g. energy-loss studies has often been emphasized [HQ. 
Proton-nucleus and deuteron-nucleus reactions were also 
used to study the Cronin-effect 0, 3 . 

Light-on-heavy nucleus- nucleus collisions offer unique 
information about the underlying dynamics, not avail- 
able in symmetric proton-proton or nucleus-nucleus col- 
lision systems. The light-on- heavy asymmetry mani- 
fests itself in an asymmetric (pcsudo-)rapidity distribu- 
tion of charged particles with respect to zero rapidity 
(or pseudorapidity) as measured by BRAHMS Q and 
PHOBOS [||. The asymmetry of the yields can be quan- 
tified by introducing the ratio of pseudorapidity densi- 
ties at a given negative pseudorapidity relative to that 
at the positive pseudorapidity of the same magnitude. 
This backward/forward ratio (forward being the origi- 
nal direction of motion of the light partner) is referred 
to as pseudorapidity asymmetry The STAR Collabora- 
tion published pseudorapidity asymmetries in 200 AGcV 
dAu collisions for several identified hadron species and 
total charged hadrons in the pseudorapidity intervals 
\r)\ < 0.5 and 0.5 < |r/| < 1.0 0. Asymmetries with 
the backward/forward ratio above unity for transverse 
momenta up to k 5 GeV/c are observed for charged 
pion, proton+anti-proton, and total charged hadron pro- 
duction in both rapidity regions. Because of the impor- 
tance of asymmetric collisions, we anticipate that proton- 
lead or deuteron-lead data will be collected at the Large 
Hadron Collider (LHC). 



The theoretical relevance of asymmetric collision sys- 
tems (in particular colliding the lightest nuclei like pro- 
tons or deuterons on a heavy partner) in testing param- 
eterizations for nuclear shadowing and models for initial 
multiple scattering was recognized prior to the availabil- 
ity of 200 AGcV RHIC deuteron-gold data @, i, EJ. 
One can take advantage of the fact that positive (for- 
ward) rapidities correspond to large parton momentum 
fractions from the light partner and small momentum 
fractions from the heavy nucleus (and vice versa for neg- 
ative rapidities), and that there are more collisions suf- 
fered traversing the heavy partner. Wang made predic- 
tions for pseudorapidity asymmetries on this basis 0] . A 
subsequent calculation focused on the lowest transverse 
momenta where pQCD may be applicable . More re- 
cently, two versions of a pQCD-based model with nuclear 
modifications were used to address deuteron-gold colli- 
sions and results were compared to available data (l3.[I3|. 
The related subject of forward-backward rapidity corre- 
lations in asymmetric systems is treated in [14| u sing a 
derivative of string percolation models, and in [15| under 
the framework of color glass condensate. 

In the present study we investigate the roles of nuclear 
shadowing and multiple scattering in the generation of 
rapidity asymmetry in the backward/forward yields at 
intermediate and high transverse-momenta. Here we use 
the HIJING shadowing parameterization and the re- 
cently released Eskola-Paukkuncn-Salgado (EPS08) nu- 
clear parton distribution functions (nPDFs) [13] • While 
the former has been applied widely, the latter was not 
available at the time of the earlier studies mentioned 
above. We calculate the pseudorapidity asymmetry in 
three representative asymmetric light-on-heavy systems: 
pBe at 30.7 GeV (Fermilab), dAu at 200 AGeV (RHIC), 
and dPb at 8.8 ATeV (LHC). We concentrate on neu- 
tral pion production and compare results to experimental 
data from the E706 experiment [3, PHENIX [lj|, and 
STAR 0. 

The paper is organized as follows: in Sec. [II] we review 
the basic formalism of the pQCD-improved parton model 
with intrinsic transverse momentum and multiple scat- 
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tering as applied to proton-nucleus (pA) and deuteron- 
nucleus (dA) collisions. This section also includes the def- 
inition of the pseudorapidity asymmetry and a discussion 
of the roles of nuclear shadowing and multiple scattering 
in asymmetric collisions. We present the results of our 
calculation in Sec. IHII and conclude in Sec. HVl 



II. CALCULATIONAL FRAMEWORK 

A. Description of the model 

The invariant cross section for the production of final 
hadron h from the collision of nucleus A and nucleus B 
(A + B-^h + X), can be written, in collincarly factorized 
pQCD. as 

E h ^-^- = E/ d2fe d * r tA(b)t B (\b- f\)dx a dx b 

abed 

dz c f a / A (x a ,Q 2 ) fb/B{xb,Q 2 ) 
da{ab-+cd) D h/c (z c ,Q 2 f ) 

T? 2~"^ sS (s+t+u) , (1) 

at nz£ 

where x a and Xb are parton momentum fractions in A 
and B, respectively, and z c is the fraction of the par- 
ton momentum carried by the final-state hadron h. The 
factorization and fragmentation scales are Q and Qf, re- 
spectively. As usual, s, t, and u refer to the partonic 
Mandelstam variables, the massless parton approxima- 
tion is used, and 



t A (s) = dzp A (s,z) 



(2) 



is the Glauber thickness function of nucleus A, with the 
nuclear density distribution, p A (s, z) subject to the nor- 
malization condition 



ft d 2 s dzp A (s, z) = A 



(3) 



The quantity da(ab — > cd)/dt in eq. |T]) represents 
the perturbatively calculable partonic cross section, and 
Dh/c{ z cQ 2 f) stands for the fragmentation function of par- 
ton c to produce hadron h, evaluated at momentum frac- 
tion z c and fragmentation scale Qf. 

The collincar parton distribution functions (PDFs) can 
be generalized to include a transverse momentum degree 
of freedom, kx, as required by the uncertainty principle. 
This can be formally implemented in terms of unintc- 
grated PDFs [I(|[2l[. To avoid some of the complications 
associated with using unintegrated PDFs, it is expedient 
to parameterize phenomenologically the kr dependence 
of the parton distributions. In this light a phenomcnolog- 
ical model for proton-proton (pp), proton- nucleus (pA), 
and nucleus-nucleus (AB) collisions incorporating parton 
transverse momentum in the collinear pQCD formalism 
was developed in Ref. (22| . In order to make the present 



study reasonably self-contained, we include details rele- 
vant to the present treatment. 

In this model, the invariant cross section can be written 



Eh^-^- =Y,f d2b d * r t A (b)t B (\b- r\)dx a dx b 

^ abed 

dk Ta dk Tb dz c f a / A (x a ,k Ta ,Q 2 ) f b /B{xb,k T b, Q 2 ) 

da(ab-^cd) D h/c (z c ,Q 2 f ) , ., 
= — J —sS(s+t+u) , (4) 



dt 



where the fey-broadened parton distribution in the nu- 
clcon is written, in a simple product approximation, as 



f a /N{x, k Tl Q 2 ) — ► g(k T ) ■ f a /N(x, Q 2 ) 



(5) 



with f a /N( x i Q 2 ) denoting the standard collinear PDF in 
the nucleon. The transverse momentum distribution is 
taken to be a Gaussian, 



exp(- 



t)pp) 



ir(k 2 



(6) 



T/PP 



where (k^} pp is the two-dimensional width of the 
transverse-momentum distribution in the proton. Based 
on the then-available pion and unidentified hadron pro- 
duction data, an estimate for the model width of the 
transverse-momentum distribution of partons in the pro- 
ton ((kj*}pp) was presented in Ref. [22j. A recent sum- 
mary of the energy dependence of the model parameter 
(fcj-), related to the average transverse momentum of the 
created pair by 



{k T ) pp 



(pt) 2 ■ 

Xf J- I pair 



(7) 



can be found in Ref. [23|]. The value of (fcy) PP increases 
logarithmically with y/s, and the data are well described 
by the function 



(PT)pair 



= (1.74 ± 0.12) . log 10 (Va) + (1-23 ± 0.2) 



(8) 



Using cq. ([8]) one can estimate (kj,) at required cms ener- 
gies. It should be noted that eq. (JSJ) is a fit to the data at 
lower energies. Since no data exists at LHC energies, we 
extrapolate (fcy). We observe that: (a) Since the depen- 
dence is logarithmic, the energy step we make from RHIC 
to LHC is comparable to the one from SPS to RHIC. We 
thus do not a priori expect a radical change in the trend, 
(b) The kr effects have been shown to be appreciable at 
low pt, and since our predictions at LHC cover a wide px 
range up to hundreds of GeV/c in pr, potential uncer- 
tainties arising from the extrapolation of Eq. 8 will affect 
a small fraction of this range at the low px end. (c) The 
present study, at LHC energies, is exploratory. Results 
of pp collisions at the LHC will be helpful in determining 
the magnitude of the presently uncertain (fcy)pp. 

It is easy to appreciate the physical necessity of the 
presence of a transverse-momentum degree of freedom 
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in proton-proton and, therefore, nuclear collisions. How- 
ever, the handling of the effect of the nuclear environment 
on transverse momenta is one of the specific features of 
the given description. Most shadowing parameterizations 
include at least some of the effects of multiple scattering 
in the nuclear medium, while the HIJING parameteri- 
zation (as we discuss further in Sec. HTC]) needs to be 
augmented with modeling nuclear multiscattering. For 
this purpose, in pA collisions we use a broadening of the 
width of the transverse momentum distribution ([6]) ac- 
cording to 



(9) 



where (k T ) pp is the width already present in proton- 
proton collisions, h p A{b) is the number of effective 
nucleon-nucleon (NN) collisions as a function of nucleoli 
impact parameter b, and C is the average increase in 
width per NN collision. In nucleus- nucleus (AB) colli- 
sions, the pt distributions of both nuclei are subject to 
the broadening represented by eq. @. 

The function h p A(b) can be written in terms of the 
number of collisions suffered by the incoming proton in 
the target nucleus, VA{b) — <JNNtA{b), where ctnn is the 
inelastic nucleon-nucleon cross sections. It was found in 
Ref. [22| that only a limited number of collisions is ef- 
fective in broadening the transverse momentum distribu- 
tion. This model scenario was referred to as "saturation" , 
and an optimal description was found with the effective 
number of nucleon-nucleon collisions maximized at 4, and 
the average width increase per NN collision, C, set to 
0.35 GcV 2 /c 2 . We do not change the values of these pa- 
rameters in the present study. The resulting transverse 
momentum broadening may appear too large relative to 
what is observed in Drell-Yan data at FNAL. However, 
we focus on meson production in the present application, 
where the PHENIX experiment extracts intrinsic trans- 
verse momenta in pp collisions at RHIC energies that are 
similar in magnitude to the ones discussed herc[24|. Fur- 
ther details about this aspect of the model can be found 
in Ref. 

The intrinsic transverse momentum kx is treated phe- 
nomenologically in Ref. [22j and in the present study. 
While next-to-leading-order (NLO) calculations provide 
a more accurate description of the parton-level cross sec- 
tion, they continue to rely on the factorization theorem 
and represent the non-perturbative information in terms 
of PDFs and fragmentation functions. Since these func- 
tions are fitted to the same data as in LO, the expected 
change is a shift in responsibilty between the perturba- 
tive and non-perturbative sectors in describing the data 
used to define the non-perturbative ingredients. In apli- 
cations to a larger set of data this will not eliminate the 
need for the phcnomcnological use of a transverse mo- 
mentum distribution, in particular considering the fact 
that even NLO may be far from a full perturbative ex- 
pansion. In addition, above we argue for inclusion of a 
transverse momentum degree of freedom on fundamental 
physical grounds as basic as the uncertainty principle. 



The simple Gaussian representation of this physics pro- 
vides a phenomenologically useful additional parameter 
(the width). This holds true both at the LO and NLO 
levels. Going to NLO may change the range of trans- 
verse momenta where intrinsic kx is important (see e.g. 
Refs. [iH [H, H3|)- No attempt is made in the present 
work to discuss the x, flavor (quarks and gluons), and ra- 
pidity dependence of kx- We have limited ourselves here 
to a simple effective description. This seems adequate for 
now, since kr effects are appreciable only at relatively low 
Pt, and a major focus of this work is asymmetry at high 
Pt (except at very forward rapidities where pt-s are low 
due to phase space constraints) . The intrinsic transverse 
momentum, (k T ) pp , enhances hadron production yields 
at low pt in both negative and positive pseudorapidity 
regions, thus tending to cancel out in the asymmetry ra- 
tio. Multiscattering, on the other hand, does impact ap- 
preciably on the asymmetry even at low pr- Overall, cal- 
culations within this framework have proven their value 
in the interpretation of hadron-production data [22l |28| . 
(For a comparison of the leading-order fcy-factorized ap- 
proach and next-to-leading order collinear approach see 
Ref. H|.) 



The collinear nPDFs f a /A( x ,Q 2 ) ar c expressible as 
convolutions of nucleonic parton distribution func- 
tions (PDFs) f a /N( x ,Q 2 ) an d a shadowing function 
Sa/A{ x i Q 2 ) which encodes the nuclear modifications of 
parton distributions. We use the MRST2001 PDFs [H 
for the nucleon parton distributions, and for the shadow- 
ing function we employ both the EPS08 shadowing rou- 
tine [lpand HIJING 0. (Other nPDFs, like FGS [H, 
HKN |32|. and the earlier EKS [IH, are used elsewhere 
to calculate pseudorapidity asymmetries [H|].) For the 
final hadron fragmentation we utilize the fragmentation 
functions in the AKK set [34[. The factorization scale 
is tied to the parton transverse momentum via Q = 
(2/3)pt/z c , while the fragmentation scale varies with the 
transverse momentum of the outgoing hadron, according 
to Qf = (2/3)pr, following the best-fit results obtained 
in Ref. [35j. This scale fixing is not modified from ear- 
lier applications, and is used consistently throughout the 
present calculations. To protect against divergences in 
the partonic cross sections, a cutoff regulator mass is nec- 
essary. We have tested the sensitivity of our calculations 
at different cms energies by varying the regulator mass 
between 0.5 — 2.0 GeV. For pt > 1.5 GeV/c our results 
show little sensitivity to variation of regulator mass. We 
present results for pt > 2.25 GcV/c here, and thus regu- 
lator mass effects are expected to be minimal. We obtain 
the density distribution of the deuteron from the Hulthcn 
wave function (36| (as in Ref. [37[), while a Woods-Saxon 
density distribution is used for gold and lead with pa- 
rameters from Ref. [38f . 
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B. Forward and backward nuclear modifications 

The nuclear modification factor is designed to com- 
pare, as a ratio, spectra of particles produced in nuclear 
collisions to a hypothetical scenario in which the nuclear 
collision is assumed to be a superposition of the appropri- 
ate number of nucleoli- nucleoli collisions. The ratio can 
be defined as a function of px for any produced hadron 
species h at any pseudorapidity 77: 



Rab(pt,v) 



1 Ehd^B^Plv 

' E h d 3 o* /d 8 p|„ 



bin) 



(10) 



where (Nbi n ) is the average number of binary collisions in 
the various impact-parameter bins. Nuclear effects man- 
ifest themselves in Rab(Pt,v) values greater or smaller 
than unity, representing enhancement or suppression, re- 
spectively, relative to the NN reference. 

In asymmetric collisions, hadron production at forward 
rapidities may be different from what is obtained at back- 
ward rapidities. It is thus of interest to study ratios of 
particle yields between a given pseudorapidity value and 
its negative in these collisions. The pseudorapidity asym- 
metry YAsymiPr) is defined for a hadron species h as 



Y Asym(PT) - E h 



d 3 - h 



a AB 



d 3 p 



E, 



d 3 p 



(11) 



Let us consider the (double) ratio of the forward and 
backward nuclear modification factors in dAu collisions 
for species h: 



R>t) 



n dAi 



(PT,-V) 



R h d Au(PT,V) 



E h d 3 a h dAu /d 3 p\_ v / E h d 3 a h dA Jd 3 p\ 
Ehd 3 a!l p /d 3 p\- v 



E h d 3 a£ p /fr 

P\v 



(12) 



As discussed in Ref. [13[ , since the pp rapidity distribu- 
tion is symmetric around y = 0, if the same backward 
and forward (pseudo)rapidity ranges are taken in both 



directions (i.e. 



< \v\ < \Vmax\), then the pp yields 
cancel in eq. (fT5|) and one obtains that the ratio defined 
in (fl~2]) is identical to the pseudorapidity asymmetry (fTTj) : 



R^Pt) 



R dAu(PT,-V) 
R h d Au(PT,V) 



(13) 



Eq. p3p is useful from the cxpcrimantal point of view 
to handle the systematic errors of data and extract the 
proper pseudorapidity asymmetry It also provides a con- 
nection between measured rapidity asymmetry and the 
nuclear modification factors. 



C. Nuclear shadowing and multiple scattering 

We expect that particle production in pA (dA) colli- 
sions will have different yields in the forward and back- 
ward directions. This is because the respective partons 



have different momentum fractions (shadowing differ- 
ences) and because the forward-going parton has to tra- 
verse a large amount of matter. We found in our earlier 
studies that the HIJING shadowing parameterization is 
particularly useful for the study of multiple scattering, 
as it requires an explicit treatment of kx broadening (see 
eq. ©) for a successful description of data [22j. In this 
Section we therefore study the interplay of shadowing 
and multiple scattering in more detail using the HIJING 
parameterization. 
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FIG. 1: (Color Online) Illustration of calculated yields at 
pr = 3.0 GeV/c in minimum bias neutral pion production 
from dAu collisions at 200 ^4GeV. The solid line represents 
HIJING shadowing with intrinsic kr in the proton, the dashed 
line is HIJING with intrinsic kr and multiple scattering, while 
the dot-dashed curve is obtained by turning off any transverse 
momentum. 

In Fig. [1] we display the neutral pion yield as a func- 
tion of pseudo-rapidity in dAu collisions at a fixed px, 
namely at 3.0 GeV/c. This illustrates how the nuclear 
effects we examine modify the asymmetry of the yields. 
The figure shows the calculated distribution using the 
HIJING shadowing parameterization without any intrin- 
sic kr (dot-dashed curve), using HIJING shadowing and 
including the intrinsic kx in the nuclcon (solid), and HI- 
JING shadowing with intrinsic kx and multiple scatter- 
ing (dashed). The parameters of eq. ^ are unchanged 
from our previous studies at midrapidity [H, [3f| , and we 
have chosen a transverse momentum value comparable to 
(kx)pp, where the various effects are clearly displayed. 

We checked that when shadowing and all nuclear ef- 
fects are turned off, the distribution is symmetric around 
midrapidity, as expected. With shadowing only, we find 
an asymmetry in the distribution: the yield at a fixed 
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negative pseudorapidty, say 77 = — 2 (Au side or "back- 
ward" ) is higher than at the corresponding positive pseu- 
dorapidity (77 = 2, d side, or "forward"). Thus the pseu- 
dorapidity asymmetry, YAsym-, is greater than unity in 
this case. The inclusion of the intrinsic transverse mo- 
mentum in the proton significantly increases the yields 
on both the Au side and the d side at this transverse mo- 
mentum, as it makes larger px-s accessible. We tested 
that adding intrinsic transverse momentum in the pro- 
ton in the absence of shadowing does not destroy the 
forward/backward symmetry. When intrinsic kx in the 
proton is added to the calculation using shadowing, but 
without the multiple scattering contribution, the sense of 
the asymmetry given by shadowing (backward/forward > 
1) is preserved. However, when multiple scattering is in- 
cluded, the yield shows a stronger increase in the forward 
direction. This is understandable in the present picture, 
since forward-going products originate in the partons of 
the deuteron, and have to traverse a large amount of nu- 
clear matter resulting in strong multiple scattering, while 
the backward products suffer no or little multiple scatter- 
ing. This has the effect of reversing the asymmetry: the 
yield on the e?-side is now greater than that on the Au- 
side. The calculated pseudorapidity asymmetry, YAsym-, 
turns out to be less than unity in this case. An alternative 
way to summarize the situation is to say that shadowing 
suppresses the yield more on the d side (forward) rela- 
tive to a symmetric collision, while the multiple scatter- 
ing contribution is understandably large in the forward 
direction. 

We have carried out similar studies at higher px values. 
When the intrinsic transverse momenta are small relative 
to the pt of the final hadron, kx effects become naturally 
smaller. Shadowing effects also become smaller as the 
antishadowing region of the HIJING parameterization is 
approached. Thus, at px > 15 GeV/c the influence of 
multiple scattering and intrinsic kx become negligible. 
These phenomena are most important at intermediate 
p T values, 2 GeV/c < p T < 8 GcV/c at RHIC. It is 
interesting to note that a similar transverse-momentum 
region is sensitive to nuclear effects at lower (e.g. CERN 
SPS) energies, due to the ~ log(y/s) scaling of the Cronin 
peak 0,1111. 



III. RESULTS 

To judge the success of the model in reproducing spec- 
tra, here we first present calculated spectra for neutral 
pion production at midrapidity and non-zero rapidities. 
This will help select the model choices providing the best 
agreement with the experimental information. The se- 
lected model variants are then used to calculate pseudo- 
rapidity asymmetries. 



A. Midrapidity spectra 

Figure [5] displays midrapidity RHIC ir° spectra from 
dAu collisions at \r/\ < 0.35, and for reference, the spectra 
from pp collisions at the same energy (-y/s = 200AGeV). 
The results are compared to data from the PHENIX col- 
laboration [l9[ . In the top left panel we show calcu- 
lated dAu spectra using EPS08 nPDFs with and with- 
out intrinsic kx in the proton, but without any multiple 
scattering contribution. The top right depicts the dAu 
spectra using HIJING plus kx with and without multiple 
scattering. The bottom panels contain data/theory (i.e. 
data/model) ratios. As can be seen clearly in the bottom 
left dAu data/theory panel, including intrinsic kx in the 
proton increases the calculated yield mostly at low px- 
We consider the EPS08 description with intrinsic kx in 
the proton satisfactory in the 4 GeV/c < px < 10 GeV/c 
interval. Comparing the two bottom data/theory panels 
it is easy to see that the HIJING parameterization gives 
a very similar accuracy for dAu when multiple scattering 
is included. This is also very close to the data/theory for 
pp collisions given by HIJING. 




p T [GeV/c] 



FIG. 2: (Color Online) Top: spectra for dAu it production 
at \n\ < 0.35. The left panel shows EPS08 with (dashed) 
and without (solid) intrinsic kx- The right panel is HIJING 
plus kr with (solid) and without (dashed) multiple scattering. 
Filled triangles denote the PHENIX data [ljj- The pp spectra 
from the same experiment and calculated with (dashed) and 
without (solid) intrinsic kx axe also included. Bottom: corre- 
sponding data/theory ratios for (left) dAu with and without 
kx and (right) pp with kx and HIJING plus kx with multiple 
scattering. 

Fig. [3] shows the spectra for neutral pion production 
from pBe collisions at \rj\ < 0.2. The top left panel 
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B. Spectra at non-zero rapidities 
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FIG. 3: (Color Online) Spectra for p + Be -> tt° + X at 
??| < 0.2. The left panel depicts EPS08 calculations with 
(dashed) and without (solid) intrinsic kr- The right panel is 
HIJING plus kr with (solid) and without (dashed) multiple 
scattering. Filled triangles denote the E706 data [13]. Data 
and calculations for proton-proton (pp) collisions are at 530 
GeV/c. 



shows the spectra using EPS08 nPDFs with and with- 
out intrinsic fcy. The top right depicts the spectra using 
HIJING plus fcx with and without multiple scattering. 
The agreement with the E706 data from Fermilab is quite 
good as can be seen from the lower panels which display 
the data per theory ratio. The HIJING parameterization 
with and without multiple scattering gives almost iden- 
tical data/theory ratios, i.e. multiple scattering has only 
a small effect on spectra from pBe collisions. This is rea- 
sonable in view of the fact that in the light Be nucleus 
there are few scattering centers. The effect of intrinsic kr 
is to increase the calculated yield, leading to data/theory 
ratios less than unity for all relevant pt- 

We conclude from Fig.s and [3] that it is necessary to 
include the intrinsic transverse momentum of partons in 
the proton in our model to obtain a satisfactory descrip- 
tion of available data. The EPS08 and HIJING parame- 
tcrizations differ to the extent that while HIJING calls for 
the inclusion of the broadening of the intrinsic transverse 
momentum distribution via multiple scattering, EPS08 
appears to incorporate this physics in their nPDFs. We 
thus concentrate on three model variants (EPS08 with 
proton intrinsic kx and HIJING with intrinsic kx and 
with/without multiscattering) in the remainder of this 
study. 



We now consider spectra at non-zero pseudorapiditics. 
We display results from EPS08 with proton intrinsic kr 
and HIJING with intrinsic kr and multiscattering. Fig. 0] 
shows the spectra and corresponding data/theory ratio 
for p + Be -> 7T° + X at -0.7 < r\ < -0.2 ("back- 
ward", left panel) and 0.2 < rj < 0.7 ("forward", right 
panel) using EPS08 nPDFs with k T and HIJING plus k T 
with multiscattering. The two sets (EPS08 and HIJING) 
give very similar data/theory ratios for both pseudora- 
pidity intervals. There is reasonable agreement with the 
E706 data [l8[ as is apparent from the quality of the 
data/theory ratios. 
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FIG. 4: (Color Online) Spectra for p+Be -> ir°+X at -0.7 < 
rj < -0.2 (left panel) and 0.2 < rj < 0.7 (right panel). The 
solid lines represent the EPS08 nPDFs with intrinsic kr, while 
the dashed line is obtained from HIJING plus kr with the 
inclusion of multiscattering. Stars denote the E706 data [l8l ]. 

The dAu spectra are displayed in Fig.[5]for both EPS08 
and HIJING at -1.0 < 77 < -0.5 (left panel) and 
0.5 < 77 < 1.0 (right panel). Note that the experimen- 
tal data from the STAR collaboration are not sep- 
arated into negative and positive pseudorapiditics, but 
rather averaged over both intervals. Therefore, the data 
points in Fig. O only serve to guide the eye; the relevance 
of the Figure is to highlight the difference between the 
EPS08 and HIJING results visible in the bottom pan- 
els. It can be seen that the EPS08 results do not dif- 
fer much in the forward and backward directions, while 
HIJING gives significantly larger data/theory ratios for- 
ward. Thus, EPS08 and HIJING differ appreciably at 
forward pseudorapiditics, 0.5 < rj < 1.0. This is not un- 
expected, because multiple scattering influences the d- 
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side (forward) more than the Au-side (backward). To be 
able to draw stronger conclusions, separated forward and 
backward data will be necessary, hopefully forthcoming 
from the high-statistics Run 8. 
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FIG. 5: (Color Online) Spectra for d + Au -> tt + X at 
— 1.0 < n < —0.5 (left panel) and 0.5 < rj < 1.0 (right panel). 
The solid line represents the EPS08 with intrinsic kr, while 
the dashed line is obtained from HIJING plus kr with the in- 
clusion of multiscattering. Stars denote STAR data averaged 
over 0.5 < \q\ < 1.0 
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FIG. 6: (Color Online) Pseudorapidity asymmetry, YAsym for 
p + Be -> 7T° + X at 0.2 < \r)\ < 0.7 (top) and 1.0 < \r)\ < 1.5 
(bottom). The solid line represents the EPS08 nPDFs, while 
the dashed line is obtained from HIJING with the inclusion of 
multiscattering. The dot-dashed line corresponds to HIJING 
without multiscattering, and filled triangles denote the E706 
data [H. 



to what is seen at lower 77. This trend is similar to what 
will be seen at other energies. 



C. Pseudorapidity asymmetry 

1. Asymmetry in pBe collisions at 30.7 GeV 

The pseudorapidity asymmetry for the rapidity inter- 
val 0.2 < \rj\ < 0.7 is shown in the upper panel of Fig. [5] 
compared with the E706 data [Hj|. The solid line rep- 
resents EPS08 with proton intrinsic &t, the dot-dashed 
line displays HIJING with intrinsic kx, and the dashed is 
HIJING plus intrinsic kr and multiple scattering. Both 
EPS08 with k T and HIJING without multiple scatter- 
ing give very small asymmetries and the data are also 
consistent with YAsym = 1 for low pr- The HIJING pa- 
rameterization with multiple scattering yields somewhat 
lower values at all transverse momenta. This is in line 
with the effect of multiple scattering moderately increas- 
ing the yield on the p-side relative to that of the i?e-sidc. 
In view of the rather large error bars, all three sets are in 
reasonable agreement with the data. The lower panel is 
our prediction for the interval 1.0 < |??| < 1.5. Here, the 
calculated effects are larger, but have a similar structure 



2. Asymmetry in dAu collisions at 200 AGeV 

Figure [7] shows the pseudorapidity asymmetry for 7r° 
production from dAu collisions at RHIC, for different 
pseudorapidity intervals. The two uppermost panels 
are our results for the asymmetry at \rj\ < 0.5 and 
0.5 < \t]\ < 1.0 compared with the STAR data 01- For 
Pt > 4.0 GeV/c, the agreement with data is quite good 
for all three sets. At lower px, multiple scattering in- 
creases the calculated yield mostly in the forward direc- 
tion as discussed in Sec. Ill CI leading to asymmetries be- 
low unity. At very high pt wc observe a divergence in 
the model predictions. The lower four panels are our pre- 
dictions for the asymmetry as pseudorapidity increases. 
The first three of these correspond to the BRAHMS pseu- 
dorapidity intervals [4l[. The general trend is that the 
asymmetry becomes larger as rj increases. This mainly 
arises from the strong shadowing in the larger nucleus 
at lower x values. Also, since increasing rj leads to de- 
creasing accessible pr due to phase space constraints, the 
effects of multiple scattering become more pronounced. 
In fact, for the largest -q considered, 3.7 < |7y| < 4.3, there 
is a marked difference between HIJING with and without 



8 



i 

0.9 
1.6 
1.4 

y 1-2 
Asym l 

0.8 

0.6 

1.6' 

1.4 

1.2 

1 

0.8 

0.6 









1 1 

111 < 0-5 - 


- /T' 

i i 4- 


i 




i 



10 



0.8<|t||<1.2 



_ I ' 


1 i 1 i 1 i 1 i . 
2.9 < |r|| < 3.5 /- 


' \ 




" i , i , i , i , i , i , r 




3 4 5 6 7 8 9 "'"2 2.5 

p T [GeV/c] 

FIG. 7: (Color Online) Pseudorapidity asymmetry, YAsym for 
d + Au — > 7T + X at different pseudorapidity intervals. The 
solid line represents the EPS08 nPDFs, while the dashed line 
is obtained using HIJING shadowing with the inclusion of 
multiscattering. The dot-dashed line corresponds to HIJING 
without multiscattering, and filled triangles denote the STAR 
data 0]. 



multiple scattering for all pt considered. 

It is pertinent at this point to make some remarks 
about pseudorapidity asymmetry in the BRAHMS pseu- 
dorapidity intervals. The BRAHMS Collaboration 
has observed a progressive suppression of both minimum 
bias nuclear modification RdAu and central-to-periphcral 
ratios, Rcp, with increasing rj. The present study is 
limited to minimum bias pseudorapidity asymmetry, and 
both shadowing paramcterizations (EPS08 and HIJING) 
adequately describe the existing experimental data at 
very forward rapidities. The EPS08 parameterization 
incorporates RHIC data at large rapidities, i.e. at low 
x, and thus reproduces the data. Therefore, at least in 
the minimum bias case, shadowing seems sufficient for 
a good description of the suppression observed at low 
x. The situation is different for the geometry-dependent 
Rep-, where (^-independent) shado wing plus conjectured 
impact parameter dependencies fLU |42| are clearly inad- 
equate in describing the observed suppression. 



3. Asymmetry in dPb collisions at 8.8 ATeV 



\rj\ < 0.9 and the lower panel is for 2.4 < \r]\ < 4.0. These 
intervals correspond to acceptance in the central detector 
and in the muon arm, respectively, of the ALICE exper- 
iment [44l ]. All three sets predict minimal asymmetry of 
the order of a few percent for the interval \r/\ < 0.9. As we 



Y 

Asym 1 
0.8 
0.6 



Y 

Asym 1 





' "*", ^ ■ 


hi < 0.9 ' 




/ 

1 . 1 


1 


I 


10 


100 300 




1 1 

s~S - — 


1 1 

2.4 < hi < 4.0 





100 300 



p T [GeV/c] 



FIG. 8: (Color Online) Predicted pseudorapidity asymmetry, 
Y A s ym for d+Pb -> tt°+X at s/s = 8.8 ATeV for |r?| < 0.9 and 
2.4 < \r/\ < 4.0. The solid line represents the EPS08 nPDFs, 
while the dashed line is obtained from HIJING with the in- 
clusion of multiscattering. The dot-dashed line corresponds 
to HIJING without multiscattering. 

move to higher 77, the predicted asymmetry becomes more 
significant. As can be seen in the lower panel of Fig. [51 
both EPS08 and HIJING predict substantial asymme- 
try up to ~ 10 GeV/c, and both variants of HIJING 
asymmetries remaining significant up to ~ 100 GeV/c, 
in contrast to EPS08. 

At the present level, neither model variant gives agree- 
ment with all aspects of the data: in an earlier calcu- 
lation we have found that shadowing paramcterizations 
which do not need to be augmented by a multiple scat- 
tering prescription [3l|, [32|, [331 have difficulty describing 
central-to-peripheral ratios at forward rapidity [12]. We 
have checked that this also holds for the EPS08 nPDFs. 
On the other hand, the HIJING parameterization with 
multiscattering yields pseudorapidity asymmetries below 
unity at low transverse momenta. This deficiency may be 
cured by allowing an ry-dependent multiscattering [45(. 



Let us now turn to our predictions for the pseudora- 
pidity asymmetry in a potential future dPb collision at 
LHC energy of 8.8 ATeV. The calculated results are dis- 
played in Fig.[SJ where the upper panel is for the interval 



IV. CONCLUSION 

Here we give a concise summary of the results of our 
calculations. We have demonstrated the usefulness of 
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asymmetric (light-on-heavy) nuclear collisions at rela- 
tivists energies. As illustrated in Sec. Ill CI the physi- 
cal differences between forward- and backward-going pro- 
duced particles arise from the different ranges of x sam- 
pled (different shadowing) and different amount of mul- 
tiple scattering. This leads to observable pscudorapid- 
ity asymmetries at some collision energies and transverse 
momenta. We have considered the effects of nuclear shad- 
owing and multiple scattering on pseudorapidity asym- 
metry for three asymmetric systems: pBe, dAu, and dPb 
in a wide energy range. 

To calibrate and fine-tune our model we first exam- 
ined spectra of produced neutral pions. We found that 
there are two avenues in the model for the reasonable de- 
scription of these data: (i) HIJ1NG shadowing, intrinsic 
kx, plus multiscattering, or (ii) some other nPDFs (we 
use EPS08 here) and intrinsic kr (but no additional mul- 
tiscattering). We then calculated pseudorapidity asym- 
metries with these prescriptions. Overall, the calculated 
asymmetries are in reasonable agreement with available 
experimental data. Intrinsic transverse momentum in the 
nucleon is seen to be important at low pr- Multiple 
scattering increases the yield in the "forward" or posi- 
tive pseudorapidity region, thus leading to a tendency for 
asymmetries less than unity at low pt in the scheme ex- 



plicitely relying on multiple scattering, at variance with 
the data. An LHC measurement at a high pseudorapidity 
and high pt (where kr effects no longer make a differ- 
ence) may be able to distinguish between strong shadow- 
ing (as in the HIJrNG prescription) and an nPDF with 
a relatively weaker gluon suppression (like e.g. EPS08). 

A major constraint in assessing pseudorapidity asym- 
metries is the limited availability of data for direct com- 
parison with theoretical calculations. More data in asym- 
metric light-on-heavy collisions separated with respect 
to positive and negative pscudorapidities are needed to 
judge calculated pseudorapidity asymmetries. At RHIC, 
it is expected that the high-statistic dAu Run 8 is going 
to provide such a large data set. 
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